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1 Introduction 



Fractional derivatives [1-6], have played a significant role in engineering, 
science, and pure and applied mathematics during the last decade [7-10]. 
Fractional calculus found many applications in recent studies of scaling phe- 
nomena [11-13]. One can find other important and interesting application of 
fractional calculus in classical mechanics [14-20]. Riewe has used the frac- 
tional calculus to develop a formalism which can be used for both conservative 
and non conservative systems [14, 15]. Although many laws of nature can be 
obtained using certain functionals and the theory of calculus of variations, 
not all laws can be obtained by this way. For example, almost all systems 
contain internal damping, yet traditional energy based approach cannot be 
used to obtain equations describing the behavior of a non conservative system 
[14, 15]. By this approach, one can obtain the Lagrangian and the Hamilto- 
nian equations of motion for the nonconservative systems . 
Klimek have introduced the symmetric fractional derivative and the Euler- 
Lagrange equations for models depending on sequential derivatives were ob- 
tained by using the minimal action principle [9]. In addition, the correspond- 
ing Hamiltonian was introduced and the Hamilton's equations were ana- 
lyzed [9] . The models described by fractional order derivatives of Riemann- 
Liouville type in sequential form were analyzed in both Lagrangian and 
Hamiltonian formalisms in [10]. 

Taking into account the constrained systems play an important role in gauge 
theories and the distinguishing feature of singular Lagrangian theories in 
the Hamiltonian formulation is the presence of constraints, we believe that 
the application of fractional calculus to constrained systems presents interest 
both from theoretical and practical point of views. Recently, an extension of 
the simplest fractional problem and the fractional variational problem of La- 
grange was obtained [18,19]. A natural generalization of Agrawal's approach 
[18,19], is to apply the fractional calculus to constrained systems [21-22] and 
to obtain the Hamiltonian formulation of constrained systems. 

The main aim of this paper is to obtain the Hamiltonian equations of 
motion for fractional variational problems with constraints. 

The plan of this paper is as follows: 
In Sect. 2 Riemnann-Liouville fractional derivatives are briefly reviewed. 
Sect. 3 is dealing with fractional Hamiltonian equations initiated by Riewe. 
Sect. 4 presents the Agrawal's approach of the fractional Euler-Lagrange 
equations. In Sect. 5 the extension of Riewe's formulation for the constrained 
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systems is obtained. Three examples of fractional constrained systems are 
investigated in Sect. 6 . Finally, Sect. 7 is dedicated to our conclusions. 

2 Riemann-Liouville fractional derivative 

Several definitions of a fractional derivative has been proposed. These defi- 
nitions include Riemann-Liouville, Grunwald-Letnikov, Weyl, Caputo, Mar- 
chaud and Riesz fractional derivatives [1-6] . Here, we formulate the problem 
in terms of the left and right Riemann-Liouville fractional derivatives, which 

are defined as [5] 

The left Riemann — Liouville fractional derivative 

1 / ^ \ " * 

- r(;r3^ U) VM*. (1) 

and 

The right Riemann — Liouville fractional derivative 

= r(^ (-5?) /(r).x, (2) 

where the order a fulfills n — 1 < a < n and F represents the Euler's gamma 
function. If a is an integer, these derivatives are defined in the usual sense, 
i-e., 

WW=(-;^)°.«=1.2..... (3) 

3 Riewe's formulation 

Now we shall briefly review Riewe's formulation of fractional generalization 
of Lagrangian and Hamiltonian equations of motion. 
Let us consider the action function of the form [14,15] 

f L{{ql,Ql,},t)dt, (4) 

J a 

where the generahzed coordinates are defined as follows: 

{a^tTXrit), Ql. mr'xrit), (5) 
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and r = 1,2,...,R denotes the number of fundamental coordinates, n — 
0, A^, the sequential order of the derivatives defining the generalized co- 
ordinates q, and n' = 1,...,N' the sequential order of the derivatives in 
definition of the coordinates Q. A necessary condition for S to admit an 
extremum for given functions Xr{t) is that Xr{t) satisfy the Euler- Lagrange 
equations [14,15] 

Following references [14,15], the generalized momenta take the form 

Pn- 2^ U^6 ) r. r ' 

k=n+l ^^k 

N' ar 

<'= E UD^'^"'-'^- (7) 

k=n'+l "^k 

The canonical Hamiltonian has the following form: 

R N-l R N'-l 

^ = E E pWu+i + E E K'QWi - L. (8) 

r=l n=0 r=l n'=0 

The Hamilton's equations of motion read as 

dH dH 

For n = 1, TV, n' = 1, N' we have the following equations of motion 

„ , OH 

^ = *D,p„, ^ = „D,7r„„ (10) 
r)M rIT 

The other equations are given by 

r)]-[ f)M 

Q-F = 'in+l = a^?q:, Q^ = Qn+l = t^tQ:.', (12) 

^ - (13) 
dt~ dt' ^^^^ 

where, n — 0, N, n' — 1, N'. 
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4 Agrawal's approach 



Agrawal have obtained the Euler-Lagrange equations for fractional varia- 
tional problems [18]. Now we would like to give review for his approach. 
Consider the action function 

S[q'„....ql^]= fL{{ql,Ql,},t)dt, (14) 

J a 

subject to the independent constraints 

^m{t, ql -, Qo^ ql Qn') ^0,m<R, (15) 
where the generalized coordinates are defined as follows: 

:= UBtTxrit), Ql, := (*Df )"'a;.(t), (16) 

Then the necessary condition for the curves Qq, ■■■■,qo with the boundary 
conditions go(^) = Qoi^) = qo'^ r = 1, 2, , i?, to be an extremal of the 
functional given by equation (14) is that the functions satisfy the following 
Euler-Lagrange equations [18]: 

^ + E it^tr^ + E UD*)"'^ - 0, (17) 

Oqo n=l ^qn n'=l ^^n' 

where L is defined as [18] 
L({g;, g;,}, t, X^t)) = L({g;, g;,}, t) + A^(t)$„(t, ql ql^, ql, Q^,), (18) 
here, the multiple Xm{t) G are continuous on [a, 6]. 

5 Fractional Hamiltonian formulation 

In fractional calculus it is not a unique way to define the Hailtonian for a 
given Lagrangian, mainly because they are several definitions of fractional 

derivatives. 

To obtain the Hamilton's equations for the the fractional variational prob- 
lems, we need to re-define the lift and the right canonical momenta as follows: 

Pn- 2^ Kt^b ) or ' 

k=n+l '^Ik 

N' 

<'= E iaJ^tr'^'-'j^- (19) 
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The canonical Hamiltonian has the following form: 



R N-1 R N'-l 



^ = E E pWn+i + E E K'QWi - L. 



(20) 



r=l n=0 r=l n'=0 



Then, the modified canonical equations of motion are obtained as 



(21) 
(22) 



Here, n = 1, N, n' = 1, N' . 

The other set of equations of motion are obtained as 

K, H} = g;^, = „Dr?;, {<„ H}} = Qi^, = ,D^g;, 



(23) 



aff _ dL 
'dt ~ ~~dt' 



(24) 



Here, n — 0, ...,iV, n' = 1, AT' and the commutator {, } is the Poisson's 
bracket and it is defined as 

, dA dB dB dA dA dB dB dA 



where, n = 0, N, n' = 1, N' . 

6 Examples 

In this section, we obtain the Hamiltonian equations for an constrained frac- 
tional variational problems. All calculations were done for < a < 1. 

6.1 Example 1 

As the first example, consider the following constrained fractional variational 
problem from the optimal control theory [18]: minimize 



(25) 




(26) 



such that 



oD^xi = - 
a;i(0) = 1. 



Xi + X2, 



(27) 
(28) 
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The modified Lagrangian L is given by 

L^^[xl + xl] + l(Pi + X(P2, (29) 

where A are Lagrange multiphers and 0i,02 are the constraints and are 
defined as 

(pi = oD^xi + xi - X2 = 0, (30) 
(/,2 = a;i(0)-l = 0. (31) 

The generahzed canonical momenta read as 

pI = 1, pI = 0. (32) 

The Hamiltonian H is given by 

H = pUI + 4] - l<t^i - (33) 

The Hamiltonian equations are obtained as 

|^=P^-/ = 0, (34) 



2 - 0, (35) 



dH 
dql 

which is identically satisfied. The other equations of motion are given by 

^-xi-l = tBUpl), (36) 
^-X2 + l^0. (37) 



dxi 
OH 



0X2 

Making use of (33), equations (36) and (37) can be put in the form, 

xi + / + iD^(/) = 0, (38) 
X2-l = 0. (39) 

Equations (38) and (39) are in exact agreement with those obtained by 
using the Euler-Lagrange formulation for fractional variational problems [18]. 
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6.2 Example 2 

Let us consider the following Lagrangian 

L= i(xi + x2)'. (40) 

The Lagrangian (j4Uj) is singular because its usual Hessian matrix has rank 1. 
One possible fractional generalization of (j4(jp is given as follows 

L = ^(oD^xi + oD^X2)^ (41) 
For Lagrangian ()41|) we introduce the notion of fractional Hessian as 



■,^,i = l,2. (42) 



Evaluating the rank of ()5H) we obtain the result being 1, so the system is a 
constrained system. 

On the other hand from (|^T| we obtain immediately the form of the 
canonical momenta as 

pI = oD^a;i + oD^X2,p^ = oD^Xi + oD^xa, (43) 
so that (jlHj) implies that 



Po-Po = 0, (44) 
which represents a primary constraint. 

Using (j4H) and (j43|) the form of the canonical fractional Hamiltonian is 
given by 

H=\{plf + ,J^'^XM-Pl)- (45) 

If we denote in (jl^ oD°X2 by A we observe that canonical fractional 
Hamiltonian becomes 

H=\{plf + Kpl-Pl)- (46) 
From ()46|) we obtain the fractional Hamiltonian motions as follows 

oD^xi=p^-A,oD°X2 = A. (47) 
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oD°pi = 0,oD>^ = 0. (48) 

Using (j44p and (j48|) we conclude that the primary surface of constraints 
is preserved after considering its evolution from the fractional calculus point 
of view. 

6.3 Example 3 

Finally, let us analyze the following Lagrangian 

L=]^{xif - X2X^. (49) 
One possible generalization of (49) is given as 

L=^(oD^a;i)2-(oD^X2)x3. (50) 
The fractional Hessian corresponding to (50) is given by 

has rank one, therefore the fractional Lagrangian ()5()|1 is degenerate admitting 
two primary constraints. 

Using (jMH) we obtain the momenta as follows 

pl = ,Btxupl = Xs,pl = 0. (52) 
From (52) we identify the primary constraints as 

pI-X3 = 0,pI = 0. (53) 
The canonical Hamiltonian in this case is given by 

H=^-^+oB1x2{pI-xs). (54) 
Denoting qD1x2 by A and by using the canonical equations become 

oD^xi = p^,oD^X3 = 0,oD^X2 = A, 

oB^pI = 0,oBtpl = -X,oBtpl = 0,pl-q3 = 0. (55) 
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We observed from (55) that the consistency conditions of the primary con- 
straint — X3 = is compatible with the equations of motion. 

The evolution of the second primary constraint makes A zero, therefore 

oT>tx2 = 0. (56) 

From (56) we conclude that a new secondary constraint appears, therefore 
the initial phase-space was modified. 

7 Conclusion 

The fractional quantization problem is an open and important issue. The 
main difficulty is related to the non-locality of fractional Lagrangian. An- 
other specific problem in this field is related to the existence of various def- 
initions of fractional derivatives which lead to different forms of fractional 
Lagrangians and fractional Hamiltonians. 

In this study, the Hamiltonian equations have been obtained for uncon- 
strained and constrained variational problem, in the same manner as those 
obtained by using the formulation of Euler-Lagrange equations for variational 
problems obtained in [18]. The Hamiltonian describes a non-conservative 
system, therefore the the Poisson brackets introduced in this study leads us 
to a Hamiltonian conserved in time in the limit a — > 1 and for Hamilto- 
nian which does not explicitly depends on time. In this paper two different 
systems with constraints were investigated, the first kind corresponds to a 
fractional Lagrangian containing Lagrange multipliers. For this case the frac- 
tional Hamiltonian and its corresponding equations were obtained. 

The second kind of constrained system involves primary and secondary 
constraints. Having in mind to obtain similar results with the usual analysis 
of constrained systems in the limit process we investigated the form of the 
equations of motion and the evolution of the surface of primary constraints 
for the last two examples. It was observed that the form of the primary 
constraints and the form of the Hamiltonian equations arc in agreement with 
each other for the second example but for the third one the consistency 
conditions impose new constraints. 

In the general case the primary constraints may produce secondary con- 
straints and the process may continue until no new constraints will appear. 
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